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EXECUTIVE  SUMMARY 


Regional  spectral  discriminants  have  been  the  subject  of  several  intensive  studies 
in  recent  years  due  to  the  establishment  of  in-country  seismic  monitoring  system  in  the 
Soviet  Union.  Moreover,  with  the  advance  made  in  high  dynamic  range,  broad-band 
instruments,  these  spectral  discriminants  may  then  be  quantified.  This  study  will  analyze 
the  discriminatory  capabilities  of  these  spectral  discriminants  using  seismic  data  collected 
from  a  system  of  state-of-the-art  seismic  instruments  installed  in  western  China.  The 
proximity  of  this  station  (WMQ)  to  the  Soviet  test  sites  in  E.  Kazakh  and  the  active 
earthquake  zones  in  the  Tien  Shan  region  makes  it  an  excellent  source  of  data  for  such  a 
regional  discrimination  study. 

Three  spectral  ratio  discriminants  are  studied,  namely  Pg  and  Lg  spectral  phase 
ratioing,  Pg  and  Lg  frequency  band  ratioing,  and  three-component  spectral  ratioing.  These 
discriminants  are  evaluated  using  a  similar  data  set  analyzed  under  a  standard  windowing 
and  processing  criterion.  The  Pg  and  Lg  spectral  ratio  performs  very  well  in  discriminating 
between  earthquakes  and  explosions  between  the  frequency  range  of  3  to  12  Hz.  The  Pg  and 
Lg  frequency  band  ratioing  does  not  show  to  be  a  promising  discriminant,  but  it  provides 
insight  to  the  understanding  of  spectra  characteristics  for  earthquakes  and  explosions.  The 
three  component  analysis  performs  very  similarly  to  the  Pg  and  Lg  spectral  ratio 
discriminant,  but  the  unpredicted  low  level  of  the  ratios  between  horizontal  and  vertical 
components  of  Lg  for  explosions  warrants  further  analysis. 

Various  parameters  that  might  have  influenced  the  characteristics  of  the  spectral 
ratios  have  been  investigated,  including  depth  effect,  frequency  effect,  and  path  effect.  By 
studying  spectral  ratios  for  earthquakes  traversing  along  paths  similar  to  those  of  the 
explosions,  it  was  concluded  that  the  spectral  ratio  discriminant  remains  effective.  Although 
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it  is  discovered  that  the  depth  of  the  seismic  events  h^s  9  controlling  effect  in  the  level  of  the 
Lg  and  Pg  spectral  ratios  in  a  kilometer  scale,  it  is  unlikely  that  such  effect  is  observable  in 
the  meter  scale.  This  study  has  shown  that  in  Central  Eurasia,  using  events  in  and  around 
the  station  WMQ,  the  spectral  discriminant  appears  to  be  an  effective  one.  To  further 
establish  this  discriminant,  studies  need  to  be  performed  in  regions  with  an  abundance  of 
ground-truth  geological  information,  thus  allowing  an  in-depth  understanding  of  the 
mcx:hanism  underlying  such  discriminant. 
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1.0  INTRODUCTION 


With  the  increasing  network  capability  of  the  seismic  monitoring  system  in  the  Soviet 
Union,  regional  discrimination  has  become  a  veiy  important  issue.  High  frequency  studies  of 
regional  phases  have  demonstrated  the  significance  of  high  frequency  regional  discriminants 
and  revealed  some  very  important  observations.  The  various  studies  of  Bakun  and  Johnson 
(1908),  Murphy  and  Bennett  (1982),  Bennett  and  Murphy  (1989),  and  Taylor  (1988)  have 
found  certain  regional  spectral  ratio  discriminants  to  be  effective  under  particular  geological 
criteria  and  presented  opposing  analysis  for  such  a  discriminant  in  other  criteria.  The 
underlying  mechanism  for  such  a  discriminant  to  perform  satisfactorily  remains  unknown  and 
may  be  attributed  to  our  lack  of  understanding  towards  the  generation  of  the  Lg  phase. 
Leading  hypotheses  of  Lg  generation  include  PS  conversion  at  the  free  surface,  scattering  of 
body  wave  energy,  and  generation  of  the  non-geometrical  5*  phase  (Gutowski  et  al,  1984). 

Short  period  regional  discrimination  relies  heavily  on  regional  phases  Pn,  Pg,  and  Lg 
which  are  well  excited  at  regional  distances  of  up  to  approximately  25°.  For  small  to 
moderate  sized  seismic  events,  these  phases  are  best  observed  on  short-period  and  broad¬ 
band  seismographs.  It  is  through  the  advance  made  in  high  dynamic  range,  broad-band 
instruments  that  these  discriminants  can  be  tested.  In  particular,  this  study  will  focus  on 
investigating  the  ability  to  discriminate  between  earthquakes  and  nuclear  explosions  through 
spectral  ratio  analysis  of  these  regional  phases  in  central  Eurasia.  The  excitation  and 
propagation  parameters  of  the  regional  crustal  phases  will  be  quantified.  Geological 
assessments  of  the  regions  of  interest  will  aid  in  the  interpretation  of  regional  variations  of 
crustal  propagation  characteristics. 

To  further  exploit  the  capability  of  the  regional  seismic  phases  Pn,  Pg,  Sn,  and  Lg  in 
discrimination  problems,  it  is  important  to  understand  fully  the  excitation  and  propagation  of 
these  phases  across  a  wide  variety  of  tectonic  and  geological  regions  for  quarry  blasts. 
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earthquakes  and  nuclear  explosions.  Discrinunants  and  yield  estimation  corrections  may  not 
be  applied  uniformly  without  regional  constraints,  owing  to  the  effect  of  varying  site 
geologies  and  propagation  paths.  In  addition  to  reflecting  the  source  characteristics,  the 
discriminants  may  also  be  affected  by  regional  variations  in  attenuation  properties,  and  in 
crustal  structure  between  each  seismic  source  and  the  various  receivers.  The  path 
attenuation  effects  are  crucial  in  seismic  yield  estimation  since  the  amplitudes  are  heavily 
dependent  on  the  efficiency  of  the  propagation  of  the  high-frequency  regional  phases. 
Propagation  of  regional  phases  generated  by  seismic  events  has  been  widely  smdied  to 
measure  their  attenuation  properties  for  corrections  to  yield  estimates.  In  particular,  Lg  is 
most  commonly  used  to  investigate  various  regions  of  the  world  due  to  its  efficient 
propagation  at  regional  distances  (c.g.,  Nuttli,  1973,  1980;  Gupta  and  Burpetti,  1980; 
Hasegawa,  1985;  Campillo  et  a/.,  1985).  In  most  of  these  studies,  the  attenuation  effects 
have  been  estimated  by  assuming  some  preferred  forms  of  the  geometrical  spreading- 
dispersion  effect  and  fitting  Q(f).  In  recent  studies  by  Nuttli  (1986a,b),  very  different 
attenuation  corrections  are  required  for  estimating  the  yields  of  NTS  and  East  Kazakh 
events  from  regional  data. 

In  general,  there  is  little  agreement  among  the  various  investigators  as  to  the  value  of 
relative  amplitudes  and  frequency  contents  of  crustal  arrivals  for  discriminating  between 
earthquakes,  quarry  blasts,  and  nuclear  explosions.  The  dominant  phase  velocities  of  Lg 
have  been  found  to  depend  on  the  source  depth  (Barley,  1979),  and  Lg  excitation  is  observed 
to  decrease  with  increasing  source  depth  compared  to  P  (Noponen  and  Bumetti,  1980). 
Murphy  and  Bennett  (1982)  claim  that  at  UBO,  explosion  spectra  have  a  lower  frequency 
content  than  earthquake  spectra.  Recently,  the  implications  of  high-frequency  regional 
propagation  characteristics  of  various  crustal  phases  with  regards  to  detectability  and 
discrimination  capability  for  decoupled  nuclear  explosions  by  a  hypothetical  internal  network 
in  the  USSR  were  discussed  by  Evernden  et  al.  (1986).  Their  study  was  based  on  certain 
assumptions  of  the  frequency  dependence  of  Q  and  assumptions  of  the  source  spectral  slope. 
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and  it  may  therefore  be  limited  in  bandwidth  and  geographical  coverage. 


Digital  seismic  data  from  the  China  Digital  Seismic  Netwoik  (CDSN)  is  utilized  to 
study  the  regional  propagation  of  high-frequency  seismic  phases  in  Eurasia  at  relatively 
close-range  distances.  This  study  utilizes  high-frequency  and  broad-band  seismic  data  from 
the  station  WMQ  in  Wulumuqi,  Xianjiang,  China  belonging  to  the  CDSN.  These  instruments 
yield  high  quality  regional  phase  data  of  up  to  20  Hz,  allowing  us  to  perform  studies  in  using 
both  spectral  ratio  and  phase  ratio  methods  for  the  different  bands  of  recordings.  The 
deployment  of  the  CDSN  network  has  made  it  possible  to  study  the  propagation  of  regional 
phases  in  Eurasia  at  distances  of  less  than  15®  (Fig.  1).  The  proximity  of  these  stations  to 
the  test  sites  in  the  Soviet  Union  and  China  has  provided  data  with  excellently  recorded  Pn, 
Sn,  Lg,  and  Pg  at  regional  distances.  These  data  constitute  an  excellent  database  for 
research  towards  understanding  the  underlying  properties  governing  the  excitation  and 
propagation  of  seismic  waves  in  the  stable  and  tectonic  regions  of  Eurasia. 
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Events  used  for  P/Lg 
in  Central  Eurasia 


70  do  90  100  110  120 


Figure  1.  Map  showing  the  locations  of  the  earthquakes  (circles)  and  E.  Kazakh  explosions 
(squares)  used  in  this  study.  The  location  of  station  WMQ  is  also  indicated.  The 
two  events  indicated  by  89043b  and  86348  in  solid  circles  are  within  20°  azimuth 
from  WMQ. 
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2.0  DATA  ANALYSIS 


Seismic  event  discrimination  based  on  spectral  signatures  has  been  applied  to  regional 
data  recorded  at  stations  of  the  CDSN  during  1987.  The  data  base  consisted  of  26  explosions 
and  101  earthquakes  (Table  1)  recorded  at  the  station  WMQ  of  the  CDSN  with  paths 
traversing  the  East  Kazakh  and  Tien  Shan  regions.  The  locations  of  the  station  WMQ  and  of 
the  seismic  events  are  shown  in  Figure  1.  In  this  study,  the  regional  phases  Pn,  Pg,  and  Lg 
are  emphasized  in  the  analysis.  Data  are  retrieved  from  the  vertical  component  of  the  short- 
period  and  broad-band  channel  of  the  CDSN  station.  In  order  to  study  the  frequency 
distribution  of  the  various  phases,  we  have  performed  a  band-pass  filtering  of  the  data  (Fig. 
2a).  For  this,  we  have  used  an  infinite  impulse  response  3-pole  butterworth  filter  in  order  to 
be  able  to  band-pass  the  entire  waveform.  The  filtered,  short-period  data  indicate  that  there 
are  Fg-phase  signals  above  the  noise  level  at  a  pass-band  of  as  high  as  16-20  Hz,  whereas 
Lg-phase  signals  are  observed  up  to  the  8-16  Hz  pass-band  for  explosions,  and  there  is 
much  lower  frequency  energy  for  earthquakes. 

Pg  and  Lg  arrival  times  are  determined  and  windowed  from  2  seconds  before  the  picked 
arrival  times  for  a  total  of  10  and  50  seconds  respectively  (Figure  2b).  The  noise  is 
determined  as  20  seconds  before  the  arrival  of  each  phase.  A  10%  cosine  taper  is  applied  to 
the  beginning  and  end  of  the  segmented  waveforms.  The  windowed  data  segments  are 
Fourier  transformed  into  the  spectral  domain.  The  spectral  ratios  are  obtained  by  taking  the 
geometric  means  at  a  particular  pass-band  for  each  phase  and  then  dividing  by  the  means  of 
each  other  phase.  Only  data  with  a  signal-to-noise  power  ratio  of  over  2  for  both  Pg  and  Lg 
phases  are  used  in  the  study.  Spectral  ratios  are  obtained  for  frequency  bands  0.2-0.5  Hz, 
0.5-2  Hz,  2-4  Hz,  4-6  Hz,  6-8  Hz  for  the  broad-band  data;  for  the  short-period  data, 
spectral  ratios  are  obtained  for  frequency  bands  2-4  Hz,  4-6  Hz,  6-8  Hz,  8.-16  Hz,  16-24 
Hz.  Time-domain  measurements  have  been  performed  for  a  few  test  cases  and  have  been 
found  to  produce  a  negligible  difference  from  those  obtained  from  the  spectral  method. 
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Figure  2a.  Band  pass  filtering  of  an  explosion  in  E.  Kazakh  recorded  at  WMQ.  An  infinite  im¬ 
pulse  response  3-pole  butterworth  filter  is  used.  The  top  trace  is  the  unfiltered 
short-period  vertical  data.  The  frequency  bands  are  0.1 -0.2  Hz,  0.2-0.4  Hz,  0.4- 
0.8  Hz,  0.8- 1.6  Hz,  1. 6-3.2  Hz,  3.2-6.4  Hz,  and  6.4-10.0  Hz. 
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3.0  RESULTS 


Guided  by  the  band-pass  filtering  analysis  which  indicates  the  broad-band 
characteristics  of  the  seismic  signals,  we  are  able  to  design  the  spectral  ratio  discrimination 
study  in  the  E.  Kazakh/Tien  Shan  region  with  the  appropriate  frequency  bands  and  phases. 
Three  discriminants  will  be  tested  in  this  study.  First,  the  spectral  ratios  between  Pg  and  Lg 
phases  will  be  examined  for  each  frequency  band.  Second,  band-ratioing  between  two 
frequency  bands  will  be  performed  to  Pg  and  Lg  phases  individually.  Third,  the  ratioing 
between  the  vertical  and  horizontal  component  Pg  and  Lg  phases  will  be  studied.  There  are 
various  parameters  that  may  affect  the  analysis,  and  these  will  be  studied  in  the  phase 
ratioing  study. 

3.1  Pg  and  Lg  Spectral  Phase  Ratioing 

The  phases  that  are  used  in  the  spectral  phase-ratioing  consist  of  Pg  and  Lg.  The 
bandpasses  used  range  from  1-2  Hz  to  18-20  Hz  for  the  short-period  data  and  from  0.5-1  to 
8-10  for  the  broad-band  data.  Examples  of  the  ratios  between  short-period  and  broad-band 
Pg  and  Lg  plotted  against  the  body  wave  magnitudes  as  obtained  from  the  NEIS  for  various 
band-passes  are  shown  in  Figures  3  and  4.  It  may  be  observed  from  the  plots  that  there  is 
distinct  separation  between  the  earthquake  and  explosion  populations  for  the  various 
passbands.  The  PgILg  ratios  for  explosions  are  over  an  order  of  magnitude  larger  than  those 
for  earthquakes.  Figure  4  shows  that  there  is  no  distinct  difference  in  spectral  ratios  between 
two  groups  of  earthquakes  that  traverse  paths  that  are  of  different  azimuths  from  WMQ. 

Spectral  ratios  between  Pg  and  Lg  plotted  against  the  frequency  for  the  short-period  and 
broad-band  data  are  shown  in  Figures  5  and  6.  For  short-period  data,  it  is  observed  that 
within  the  frequency  range  of  2  to  14  Hz,  the  spectral  ratios  between  Pg  and  Lg  phases  are 
capable  of  discriminating  between  earthquakes  and  explosions.  In  the  case  of  broad-band 
data,  there  is  separation  between  earthquakes  and  explosions  from  4  to  10  Hz.  A  histogram 
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Figure  3a.  Spectral  ratios  of  PgILg  vs.  mb  for  2.0-4.0  Hz  using  short-period  vertical  compo¬ 
nent  WMQ  data.  The  earthquake  and  explosion  populations  are  well  separated. 
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Figure  3b.  Spectral  ratios  of  Pg/Lg  vs.  mb  for  4.0-6.0  Hz  using  short-period  vertical  compo¬ 
nent  WMQ  data.  The  earthquake  and  explosion  populations  are  well  separated. 
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Figure  4a.  Spectral  ratios  of  PgiLg  vs.  mb  for  3.0-5.0  Hz  using  broad-band  vertical  compo¬ 
nent  WMQ  data.  The  earthquakes  are  indicated  as  west  of  WMQ  or  cast  of 
WMQ.  The  earthquake  and  explosion  populations  are  well  separated. 
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Figure  4b.  Spectral  ratios  of  Pg/Lg  vs.  mb  for  4.0-6.0  Hz  using  broad-band  vertical  compo¬ 
nent  WMQ  data.  The  earthquakes  are  indicated  as  west  of  WMQ  or  east  of 
WMQ.  The  earthquake  and  explosion  populations  are  well  separated. 
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Figure  5.  Plot  of  Pg/Lg  spectral  ratios  versus  frequency  for  earthquakes  and  explosions  using 
the  vertical  component  short-period  data  of  WMQ.  The  two  events  within  20®  azi¬ 
muth  from  WMQ  are  indicated  as  pluses  and  diamonds. 
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Figure  6.  Plot  of  Pg/Lg  spectral  ratios  versus  frequency  for  earthquakes  and  explosions  using 
the  vertical  component  broad-band  data  of  WMQ. 
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plotted  in  Figure  7  shows  that  at  a  center  frequency  of  7  Hz,  the  earthquake  and  explosion 
populations  are  distinctly  separable  from  each  other. 

3.2  Pg  and  Lg  Frequency  Band-Ratioing 

It  has  also  been  noticed  that  earthquakes  and  explosions  have  different  energy  content 
at  different  frequency  bands.  These  variations  in  frequency  bands  directly  contribute  to  the 
possible  discriminant  using  regional  phases.  In  order  to  further  study  the  frequency  content  of 
the  Pg  and  Lg  phases,  the  spectral  ratios  at  various  frequency  bands  are  studied  to  examine 
the  spectral  characteristics.  Due  to  the  broad-band  nature  of  the  CDSN  stations,  we  are 
able  to  examine  spectral  ratios  between  widely  spaced  frequency  bands  between  1. 6-2.5  Hz 
and  16.0-25.0  Hz  using  the  short-period  data. 

Two  different  frequency  bands  of  Pg  phase,  2-4  Hz  and  4-6  Hz,  were  ratioed  using  the 
short  period  data  and  the  results  are  shown  in  Figures  8a  and  8b.  The  earthquake  and 
explosion  populations  are  mixed  with  each  other  and  show  the  same  increasing  trend  with 
magnitude.  A  similar  feature  is  observed  for  the  frequency  band-ratioing  between  2-4  Hz 
and  8-10  Hz  for  the  Pg  phase  using  short-period  data.  The  results  here  indicate  that  for 
larger  explosions  or  earthquakes,  the  Pg  phase  appears  to  have  less  high-frequency  energy. 
The  rate  of  deficiency  of  high  frequency  energy  with  respect  to  magnitude  remains  the  same 
for  the  frequency  bands  4-6  Hz  and  8-10  Hz. 

For  Lg  phase,  two  similar  bands  of  ratioing  were  performed  and  the  results  are  shown  in 
Figures  9a  and  9b.  For  the  frequency  band-ratioing  between  2-4  Hz,  a  pattern  similar  to 
that  of  the  Pg  phase  is  observed  indicating  that  there  is  a  gradual  decrease  of  high  frequency 
energy  in  the  band  of  4-6  Hz  compared  to  those  at  2-4  Hz  as  the  body  wave  magnitude 
increases.  Frequency  band-ratioing  between  2-4  Hz  and  8-10  Hz  show  a  steeper  slope 
than  the  ratios  studied  earlier,  shown  in  Figure  9a.  This  result  indicates  that  for  Lg,  there  is  a 
more  rapid  depletion  of  energy  with  body  wave  magnitude  as  the  frequency  increases. 
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Figure  7.  Histograms  showing  the  distribution  of  PILg  log  ratio  for  earthquakes  and  explo- 
sions(shaded)  at  a  center  frequency  of  7  Hz. 
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Figure  8a.  Plot  of  Pg(2-4  Hz)//*g(4-6  Hz)  spectral  ratios  versus  mb  fw  earthquakes  and 
plosions  using  the  vertical  component  short-period  data  of  WMQ. 
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Figure  8b.  Plot  of  Pg(2-4  Hz)/Pg(8-10  Hz)  spectral  ratios  versus  mb  for  earthquakes  and 
explosions  using  the  vertical  component  short-period  data  of  WMQ. 
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Figure  9a.  Plot  of  Lg(2-4  Hz)/Lg(4-8  Hz)  spectral  ratios  versus  mb  for  earthquakes  and  ex 
plosions  using  the  vertical  component  short-period  data  of  WMQ. 


22 


Figure  9b.PIot  of  Lg(2-4  Hz)/Lg(8-10  Hz)  spectral  ratios  versus  mb  for  earthquakes  and  ex 
plosions  using  the  vertical  component  short-period  data  of  WMQ. 


3.3  Three-component  Analysis 


The  CDSN  stations  are  equipped  with  3-component  seismometers.  The  ratioing 
between  the  horizontal  and  vertical  component  data  may  guide  us  in  determining  the 
characteristics  of  the  3-comp6nent  data  and  to  evaluate  whether  the  3-component  data  may 
prove  to  be  an  effective  discriminant.  The  horizontal  data  are  first  rotated  to  be  radially  and 
transversely  polarized.  The  vertical  component  energy  represents  the  Rayleigh-type  motion, 
whereas  the  energy  on  the  transverse  component  represents  that  of  the  Love-type  motion. 

The  frequency  phase  ratioing  between  the  transverse  and  vertical  component  of  the  data 
indicates  that  the  ratios  for  the  explosions  are  less  than  those  of  the  earthquakes  and  less 
than  unity,  in  particularly  at  frequency  bands  of  1. 0-2.0  and  2.0-4.0  Hz,  as  is  shown  in 
Figures  10a  and  10b.  This  observation  is  opposite  from  what  one  may  expect  for  explosions. 
It  is  commonly  believed  that  explosions  possess  less  Love-type  energy  and  therefore  the 
ratios  Lg{Z)ILg(J)  should  be  greater  than  1.  The  ratios  between  horizontal  and  vertical 
components  of  Lg  may  not  be  used  as  any  kind  of  discriminant,  but  the  low  level  of  the  ratios 
for  explosions  is  rather  intriguing. 

We  have  further  tested  the  3-component  frequency  phase-ratio  between  Pg  on  the 
vertical  component  and  Lg  on  the  transverse  component  at  various  frequency  bands.  In 
Figures  11a  and  11b,  the  explosion  and  earthquake  phase  ratios  are  plotted  versus  body- 
wave  magnitude,  mb.  Similar  to  the  analysis  in  the  previous  section,  the  Pg{Z)  to  Lg{T) 
ratios  for  the  earthquakes  are  lower  than  those  of  the  explosions  for  both  frequency  bands 
between  6.0-8.0  and  4.0-6.0  Hz.  The  3-component  phase  ratioing  has  produced  the  same 
result  as  the  vertical  component  ratioing. 
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Figure  10a.  Plot  of  Lg(Z)/Lg(T)  spectral  ratios  at  1.0-2.0  Hz  versus  mb  for  earthquakes  and 
explosions  using  the  broad-band  data  of  WMQ. 
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Figure  10b.  Plot  of  Lg{Z)ILgiT)  spectral  ratios  at  2.0-4.0  Hz  versus  mb  for  earthquakes  and 
explosions  using  the  broad-band  data  of  WMQ. 
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Figure  11a.  Plot  of  PgiZ)ILg(J)  spectral  ratios  at  6.0-8.0  Hz  versus  mb  for  earthquakes  and 
explosions  using  the  short-period  data  of  WMQ. 
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4.0  DISCUSSION 


The  analysis  in  this  study  has  demonstrated  the  characteristics  of  spectral  ratios 
between  Pg  and  Lg  phases  for  paths  traversing  Central  Eurasia.  It  is  observed  that  the 
results  indicated  the  possibility  of  using  Pg  and  Lg  spectral  ratios  as  a  discriminant  between 
earthquakes  and  explosions.  Before  one  can  be  conclusive  on  such  a  discriminant,  a  suite  of 
parameters  need  to  be  evaluated  in  terms  of  how  they  might  affect  the  discriminant. 

4.1  Depth  Effect 

For  most  of  this  study,  the  depths  of  the  earthquakes  used  have  been  limited  to  less 
than  28  km.  This  limit  is  arbitrarily  chosen  to  ensure  that  all  earthquakes  be  of  shallow  focus. 
The  gross  effect  of  depths  on  the  spectral  ratio  relationship  is  shown  in  Figures  12a  and  12b. 
The  figures  indicate  that  the  deeper  focus  earthquakes  have  a  larger  PgILg  spectral  ratio  at 
both  the  low  (0.2-0.5  Hz)  and  the  high  (2.0-4.0  Hz)  frequency  bands.  This  analysis  indeed 
shows  that  the  depth  parameters  do  indeed  have  significant  effects  on  the  PgILg  spectral 
ratio  level.  Furthermore,  by  plotting  the  spectral  ratios  against  depth  as  is  shown  in  Figure 
13,  it  is  observed  that  the  explosion  population  remains  well  separated  from  the  earthquakes 
at  all  depths.  Frankel  (1989)  had  demonstrated  the  effects  of  source  depth  on  the  spectral 
characteristics  of  regional  phases  using  synthetic  seismograms.  Even  though  the  depth 
wavelength  in  Frankel’s  study  is  orders  of  magnitude  different  from  the  above  analysis,  it 
nonetheless  demorisirates  the  large  scale  effect  of  source  depth  on  spectra. 

4.2  Frequency  Effect 

It  is  also  a  goal  of  this  study  to  understand  the  effect  of  frequency  on  the  spectral  ratio 
discriminant.  Section  3.2  showed  the  effect  of  frequency  band-ratioing  of  the  Pg  and  Lg 
phases,  indicating  that  the  variations  in  frequency  bands  are  too  insensitive  to  allow  the 
band-ratioing  method  to  be  an  effective  discriminant.  Due  to  the  scatter  in  the  band-ratioing 
for  both  phases  (Figures  8  and  9),  it  fails  to  provide  a  definite  indication  as  to  whether 
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Figure  12a.  Plot  of  Pg/Lg  spectral  ratios  at  0.2-0.5  Hz  versus  mb  for  earthquakes  and  explo¬ 
sions  using  the  vertical  component  broad-band  data  of  WMQ. 
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Figure  12b.  Plot  of  Pg/Lg  spectral  ratios  at  2,0-4.0  Hz  versus  mb  for  earthquakes  and  explo¬ 
sions  using  the  vertical  component  short-period  data  of  WMQ. 
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Figure  13.  Plot  of  Pg/Lg  spectral  ratios  at  2.0-4.0  Hz  versus  depth  for  earthquakes  and 
plosions  using  the  vertical  component  short-period  data  of  WMQ. 


eartliquakes  or  explosions  contain  more  high  frequency  energy  for  these  two  phases.  One 
exception  is  that  there  is  an  indication  that  explosion  Pg  contains  more  high  frequency 
energy,  as  is  shown  in  Figure  8b  where  most  explosion  ratios  fall  below  those  of 
earthquakes.  In  the  case  of  Lg,  the  results  are  too  scattered  to  reach  any  conclusion. 

4.3  Path  Effect 

It  has  been  argued  that  path  effect  is  a  rather  important  factor  in  spectral  ratio  studies 
(e.g.  Frankel,  1989;  'I'aylor,  1988).  In  order  to  access  this  effect,  we  have  examined 
earthquakes  from  two  different  azimuths  from  WMQ  to  see  if  there  is  any  detected  difference 
in  their  spectral  ratios  characteristics.  Earthquakes  in  Figure  1  are  separated  into  two 
groups,  west  and  east  of  WMQ.  This  separation  is  based  upon  the  fact  that  the  region  west 
of  WMQ  may  be  more  tectonically  active  than  the  region  to  the  east  of  WMQ,  which  is  more 
nearly  a  stable  basin.  The  two  groups  of  earthquakes  are  plotted  in  Figures  4a  and  4b  in  two 
different  symbols.  There  is  a  slight  indication  that  the  earthquakes  to  the  west  may  be 
characterized  by  higher  ratios  than  those  to  the  east.  This  is  inconclusive  due  to  the  high 
degree  of  scatter  in  the  data. 

Another  approach  is  to  determine  whether  signals  propagating  along  the  path  from 
WMQ  to  E.  Kazakh  exhibit  characteristics  of  large  spectral  ratios.  We  have  identified  two 
earthquakes  that  are  within  20  degrees  in  azimuth  from  WMQ  to  E.  Kazakh  for  which  the 
signals  travel  along  a  similar  path  as  do  those  from  the  explosions.  These  events  are 
indicated  as  89043b  and  86348  in  Figure  1.  The  Pg/Lg  spectral  ratios  of  these  two  events  are 
computed  and  plotted  in  Figure  5  along  with  the  rest  of  the  results.  These  ratios  fall  among 
the  earthquake  population  and  are  among  the  largest  of  all  earthquakes,  yet  they  are  still 
distinctly  separated  from  the  explosion  population.  The  results  are  not  unexpected,  as  the.se 
two  events  are  of  uncertain  depth  (NEIS  depth  of  33km).  As  indicated  in  earlier  analysis, 
deeper  earthquakes  are  characterized  by  large  spectral  ratios.  If  these  two  events  are  indeed 
shallow  events,  the  spectral  ratios  of  these  two  events  may  serve  as  the  upper  bound  for  the 
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earthquake  populations.  If  on  the  other  hand  they  are  deeper  earthquakes,  this  fact  shall  be 
easily  discernible  and  it  will  not  be  necess^  for  such  a  discriminant  to  be  applied,  as  this 
discriminant  is  designed  for  distinguishing  shallow  focus  earthquakes  from  explosions, 
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5.0  CONCLUSIONS 


A  spectral  ratio  method  can  likely  be  used  to  discriminate  between  earthquakes  and 
explosions  in  the  Central  Eurasia  region,  but  the  lack  of  ground  truth  information  remains  an 
obstacle.  The  narrow  bandpass  filtering  scheme  provides  a  full  characterization  of  each 
seismic  signal  in  the  time  and  frequency  domains.  If  the  same  event  is  observed  at  several 
stations,  then  we  shall  have  some  descriptions  of  the  distance  dependence  of  the  wavetrain 
as  well.  The  data  characterizing  the  events  will  consist  of  descriptions  of  the  absolute  and 
relative  amplitudes  of  the  arrivals  in  various  frequency  bands  and  their  time  domain 
envelopes  (including  that  of  their  codas  as  well).  With  the  collected  data,  we  shall  also 
compile  information  that  may  be  useful  for  providing  physical  explanations  for  the  empirical 
results.  Such  information  will  include,  but  not  be  limited  to,  source  depths,  types  of  blasts, 
path  characterization,  and  source  mechanisms.  Additional  emphasis  will  be  placed  upon 
paths  that  may  have  anomalously  efficient  propagation  of  the  regional  phases. 
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